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We argue that long-range photon-photon attraction induced by the dipole inter- 
action of two electron-positron loops can lead to "vacuum self-focussing" of very 
intense laser beams. The focussing angle Of is found to increase with the beam in- 
£^ ■ tensity / as 6f ~ I 4//3 ; for the laser beams available at present or in the near future, 

•■s ; e F ^ io- 10 + io~ 7 . 

hi, 

43 ■ I. INTRODUCTION 

•rH , 

The response of quantum vacuum to strong external fields is a problem of great funda- 
mental interest. Recent technological developments made available very intense laser beams, 



and opened t 



review, see jl]). In particular, mu 
tally in light-by-light scattering 



re possibility to study non-linear QED interactions in experiment (for a recent 



ti-photon effects have been already identified experimen- 
2|. Theoretical studies of light-by-light scattering so far 
have been limited to the one-fermion-loop effects - the scattering of photons induced by the 
coupling to one electron-positron pair. The studies of higher loop corrections so far have 
been limited to two-loop effects arising from an additional internal virtual photon interact- 
ing with the single fermion loop; the corresponding results are available in an integral form 



Q; see [l] for a comprehensive review . 

The expansion in the number of fermion loops clearly converges very fast, since additional 
fermion loops bring into the effective low-energy interaction extra powers of a 2 F 2 /m 4 <C 
1, where F is the field strength tensor, m is the electron mass and a ~ 1/137 is the 
electromagnetic coupling. We thus come to the conclusion that one can safely approximate 
the photon-photon scattering amplitude by a one-fermion-loop result. 

This conclusion however may be premature for the field of an intense laser. Indeed, let us 
compare the photon-photon scattering amplitudes at the one-fermion-loop level (Fig. Ufa)) 
and two-fermion-loop level (Fig. [U(b)). The amplitude of the box diagram of Fig. U^a) in 
coordinate space falls off exponentially at large separations r as ~ (a 2 /m 4 r 7 ) exp(— 2mr) 
since it involves the exchange of a massive fermion pair in the t— channel. On the other 
hand, the amplitude of Fig. [H(b) behaves in coordinate space like ~ a 4 /m 8 r 11 (we will show 
this by an explicit calculation below). The comparison of the amplitudes of Fig. Ufa) and 
Fig. \H[b) thus leads us to the conclusion that the two-fermion-loop amplitude of Fig. [T^b) 
will dominate over the one-fermion-loop amplitude of Fig. HJa) at distances larger than 

r ~ — In a/m. (1) 

One may wonder whether multi-photon interactions can also be important. There are 
two types of such interactions: i) when additional external beam photon lines are attached 
to a single fermion loop; and ii) when additional external beam photon lines are attached to 
different fermion loops connected by t— channel photon pair exchanges. The diagrams of type 
i) for n external beam photon lines are suppressed by (eF/m 2 ) n and fall off exponentially as 
a function of inter-photon separation; we thus expect that such interactions can be neglected 
when the field strength is smaller than the critical one, F 2 ~ m 4 /a. The diagrams of type ii) 
are characterized by a power fall-off and so are potentially relevant even at sub-critical field 
strength; however the interactions involving I beam photons are suppressed by (a 3 /m 4 r 4 )'. 
This parameter becomes of the order of unity at distances smaller than r ~ a 3//4 /m. In 
this paper we consider the dynamics of laser beams of smaller intensities, corresponding to 
larger inter-photon distances (note that — lna 3> o?l 4 for a = 1/137), see Eq.Q. 



Some non-perturbative effects such as those caused by the existence of bound electron-positron states 
have also been considered lal. 
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This situation is completely analogous to the interaction of atoms: at short distances, the 
interaction involving the exchange of constituent electrons dominates, but at the distances 
large compared to the size of the atoms, the interaction is dominated by the (attractive) 
two-photon exchange despite the fact that it is of higher order in the coupling a. The key 
question now is whether the conditions for the dominance of two-fermion-loop diagrams of 
Fig. [H(b) can be realized in realistic laser beams. 

Consider a laser beam which is characterized by the wavelength A, coherence length l c 
and the transverse diameter d. The photon density in such a beam is determined by the 
ratio of the field energy density to the energy of a single photon 

Let us approximate the beam by a cylinder of height I and cross sectional area 5* with the 
symmetry axis parallel to the beam direction. The beam energy contained in the cylinder 
is then E 2 Sl. The beam intensity is the flux of energy through the unit surface, i.e. 

E 2 SL 



' S{ljc) 

hence 



E 2 c, (3) 



For example, for the Omega-EP laser at the Laboratory for Laser Energetics of the University 
of Rochester the parameters are [fjj], / = 3 • 10 20 W/cm 2 at the wavelength A = 1054 nm 
(corresponding to uo = 2tcc/X = 1.8 • 10 15 s _1 ). This amounts to the mean photon density of 
p = 5.3 • 10 28 cm -3 . Therefore the typical inter-photon distance in the beam is R ~ p -1 / 3 
which is R = 2.6 ■ 10~ 3 nm for the Omega-EP laser. 

As long as the photons obey the linear Maxwell equations the laser beam is a coherent 
state. However, the coherent state can be destroyed by the mutual interaction of the photons. 
Such interactions appear as a result of quantum fluctuations of a photon into a virtual e + e~ 
pair. These pairs can be considered as electric or magnetic dipoles of size A e = h/m e c = 
3.8 ■ 1CT 4 nm which can interact with virtual dipoles created by other photons in the beam. 
Thus, there are four dimensionful scales associated with the problem of photon interaction 
in the laser beam. These are 

A e < R < A < L , (5) 
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where L is the macroscopic length of the beam (l c or d). By comparing ([5]) and ([T]) we 
conclude that (i) the condition of ([1]) for the dominance of two-fermion-loop interactions 
can be met in the available laser beams, and (ii) since the photon wavelength A is the largest 
(microscopic) scale in the problem we can treat the two-photon interaction as a low-energy 
scattering problem. 




(a) (b) 

FIG. 1: The amplitudes of light -by-light scattering at the one-fermion loop (a) and two-fermion- 
loop (b) level. 



II. PHOTON-PHOTON INTERACTION AT LARGE DISTANCES 

Let us now turn to the calculation of the photon-photon potential density which can be 
expressed through the scattering amplitude A4(s,Q 2 ) as 

V(v,t) = - J 0^e~ iQ - x KmM(s,Q 2 ), (6) 

where x = (t, r) is the relative four-coordinate. This formula gives the potential density 
as a function of time. If we were considering the scattering of massive particles, we could 
have chosen their center-of-mass rest frame in which the potential would not depend on time 
t. This is not possible in the case of photon-photon interaction. Instead, we can obtain a 
meaningful quantity by averaging the potential over the period of the electromagnetic 
wave oscillation T = 2ix/oj. We thus define the mean potential density as 

V{T) = ±J T V{T,t)dt. (7) 

To calculate the two-photon interaction amplitude, Fig. QJb), we use the t-channel dis- 
persion relation 

1 f rIM 2 

M(s, Q) = - J M2 _Q 2 _ iQ ImX(s, M 2 ) , (8) 
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where M is the invariant mass of two t-channel photons. The imaginary part of the amplitude 
can be calculated by making the cut shown by the dotted line in Fig. [TJ(b) . This amounts 
to placing the cut photons on mass shell. According to the Mandelstam-Cutkosky rule the 
cut photon's propagator is replaced by 



2 + ;o -2^(P 2 ) £ e; V . (9) 

polariz. 

Therefore, in the low energy limit 

lmM(0,Q 2 ) = J ^6( q y((q-Q) 2 )J^\M h + ^AM 2 )\\ (10) 

where we used the fact that all photons in the beam have the same polarization (assigned 
+ here). 

The low energy on-shell amplitudes corresponding to the box diagram Fig. [TJ(a) can be 
found elsewhere (see e.g. in |7|): 

Taking the phase space integral in (TlOl) and substituting it together with (TTTj) into ([8]) yields 
the amplitude in the t-channel. However, we are interested in the s-channel amplitude. It 
can be obtained by analytic continuation to the region of space-like momentum Q. Therefore, 
the s-channel amplitude is 

, ^ 157 / a 2 V f dM 2 o 

where now Q 2 < 0. The potential density can be calculated by substituting ffl2|) into (Q. 
We have 

Integration over Q yields 



/ 



1 , o. M 2 



d 4 Q e iQ ' x — — = 4vr 2 z = K 1 (J(-x 2 + iO)M 2 ) , (14) 



where K\ is the modified Bessel function. Now we have to average the potential density (TT5|) 
over the oscillation period according to ([7]). Since the field oscillation period T = A is much 
larger than the inter-photon distance R (the typical exchanged mass is M ~ R^ 1 ), we can 
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set the integration limits in ([Tj) to infinity. The integral over t can be done by rotating the 
integration contour by tt in the plane of complex t. The result is 

V(r) = 1 151 (j*)' r iM H;*^e--<*. (15) 
1 ; (2vr) 4 8tt 2 \45m 4 j j 2T r V ; 

After integration over M we finally derive 

u 157 / a 2 \ 2 2tt 2 
™ = "^uS? (iS?) ^ 2 ' 9! (16) 

This is the mean potential density between two photons at large distances. 

Let us compare this to the corresponding result for the interaction between atoms, where 

the potential falls off as ~ r~ 6 at distances larger than the size of the atom a but smaller 

than a/ a ("Van der Waals interaction") and as ~ r~ 7 at distances larger than a/ a due 

to relativistic retardation effects ( "Casimir-Polder interaction") 2 . Our result falls off even 

faster, as r -11 , which can be traced back through the calculation to the gauge invariance of 

QED and the fact that the dipoles in our case are virtual, unlike in the interaction of atoms. 



III. DISPERSION RELATION OF THE BEAM AND SELF FOCUSSING 

In this section we would like to calculate the refraction index n of the intense laser beam 
taking into account the photon-photon attraction as described in the previous section. As- 
suming the dominance of two-particle interactions, the Hamiltonian describing the dynamics 
of interacting particles can be represented in the secondary quantization approach in the 



following general form (see e.g. 10||): 



~ Yl u p ^p^p + 2 U a p a p a P a p • ( 17 ) 

p 

The mean effective potential U for scattering of two photons corresponding to the density 
f|T6|) is given by: 

U = f d*rV(r) = ^ fjt J d 3 r J ^ e"^ M(0,Q>) 

u 157 / a 2 \ 2 f 1/R2 _ uo 157 / a 2 V 1 
= n\ r / dM 2 M s = — -. (18) 



2 The dispersive method which we use is analogous to the one used previously to describe the interaction 
of atoms and dipoles in RefsQ, P 
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The same result can be obtained directly by taking the four-dimensional integral of the right- 
hand-side of (j!4p . Denoting by Arc A the forward elastic scattering amplitude we obtain from 
(HE} U= -Aw. 

By analogy with the treatment of superfluidity we consider the quasi-particle fluctuations 
d p = c p + d p and dp = cp + a' where the background classical laser field is described by 
the commuting c— number operators c p and cp. Expanding the Hamiltonian up to quadratic 
terms in d p 's and diagonalizing the Hamiltonian yields 

H = coN - ~N NcuA + ^e(u;)6 p 6 p , (19) 

p 

where 6 p 's are the quasi-particle operators and Nq is a mean number of particles in a sphere 
of radius R: Nq ~ N • R 3 /V ~ pR? ~ 1. The dispersion law is given by 

e{u) = - A N /2 ~ uy/l - A/2. (20) 

The redshift of the laser beam occurs since a fraction of the beam's kinetic energy density 
E 2 = poo (see 02])) gets transformed into the interaction energy density (NN /2) ■ U/V = 
U p/2, where N is the total number of photons in the beam. Clearly, this fraction is \U\/2u = 
A/2. 

The phase velocity of the laser beam reads 

u = w 1 . 21 

p 4 

Eq. (|T8|) implies that A ~ _R -8 ~ p 8 / 3 . The density of the beam p decreases toward the 
beam edges. Therefore the beam phase velocity u increases towards the beam periphery. 
This implies that the photon beam will self-focus: the refraction index n = u~ x is largest 
in the center of the beam which has the "focussing lens" effect. The focussing angle can be 
determined from the following simple observation. Focussing means that the photons in the 
beam center having phase velocity 1 — A/4 will arrive to the focussing point simultaneously 
with the photons at the beam periphery having unit phase velocity. This can only happen if 
the relative ("focussing") angle between these bunches of photons satisfies cos9p = 1 — A/4. 
Since the focussing angle is small we have 

, p = = (i^Y'V^Y (22) 

v 1 \16ii 3 J \lS0m 4 R V 
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For the Omega-EP laser beam we have the following estimate: 6p ~ 1 • 1CF 10 . We 
expect an increase in the focussing angle with the increase of the beam intensity according 
to the law Op oc i?~ 4 oc p 4 / 3 oc J 4//3 , see (1181) and (HI). The maximal possible effect can 
be achieved at R ~ A e in which case Op ~ 1 • 1CF 7 . This value should be considered as 
an upper bound on the self-focussing angle due to the mechanism considered in this paper 
since our approach breaks down for R < A e . We have to emphasize that these estimates 
are admittedly qualitative, since they are derived in the mean field approximation in which 
R ~ p~ l '\ 

Our calculation shows that the focussing effect is present even for the plane wave, for 
which the leading order contribution from the box diagram of Fig. 1(a) vanishes. The 
magnitude of the expected effect however makes its experimental observation challenging. 
In particular, in a realistic experimental setup it would have to be distinguished from the 
non-linear effects caused by the presence of the air. 
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